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ABSTRACT 



A full-scale model of a double-wedge, simulating a jet 
vane used in guided missile thrust vector control (TVC) , was 
constructed for subsonic wind tunnel testing to evaluate the 
feasibility of utilizing a thermal imaging (infrared) system 
vice thermocouples to determine a surface heat transfer coeffi- 
cient. The model was preheated, then allowed to cool via 
forced convection to ambient conditions. Surface temperature 
readings were taken using thermocouples and a thermal imaging 
system during the cooldown process. These readings were then 
reduced to solve for a heat transfer coefficient (h) , using 
Newton's Law of Cooling. The results were compared with the 
theoretical solution for- a flat plate and with the results 
obtained by testing an actual TVC vane configuration under 
the identical flow conditions. 
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I. 



INTRODUCTION 



A . BACKGROUND 

The use of thrust vector control (TVC) in future cruise 
missile guidance and control systems is presently being inves- 
tigated. Removing external aerodynamic steering surfaces 
(fins) would allow installation of a single weapon configuration 
in multiple platforms with little or no restructuring of each 
platform's unique launch capabilities, thereby reducing design, 
development, and construction costs. In addition, the maneu- 
vering envelopes required of present and future missiles niay 
contain regions in which external aerodynamic surfaces are 
inadequate for control purposes. 

The incorporation of TVC systems for missile steering re- 
quires the design of jet vanes, or other deflection devices, 
that are mounted in the rocket exhaust flow. Such devices can 
be rotated to various angles of attack in order to redirect the 
flow, thus maneuvering the missile in the desired direction. 

The vanes must be made to withstand severe steady and transi- 
ent thermal loads, shock wave impingement, and erosion from 
motor exhaust particulates during extended use. 

One of the first steps in the design of such a vane is to 
develop a model of the thermal environment to which the vane 
will be subjected. This includes the determination of surface 
heat transfer coefficients for particular vane configurations 
under forced convection conditions. 
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The following discussion concerns the development of 
experimental methods to determine a heat transfer coefficient 
on a heated surface in subsonic flow. Utilizing state-of-the- 
art infrared technology, in place of thermocouples, to determine 
surface temperature gradients as a body is heated (or cooled, 
as in this case) would reduce' the time and effort required to 
develop instrumented models. The method also has the potential 
to reduce or eliminate the need for elaborate data acquisition 
systems, while maintaining an accuracy that is at least as 
good as that obtained from conventional methods. Finally, 
the methods- described here provide a global as well as local 
view of surface heat transfer. Thus it is possible to detect 
and analyze regions of particular interest that might otherwise 
go unnoticed in the data produced by point temperature 
measurements . 



B. OBJECTIVES 

The objectives of this study are to: 

1) Design a typical TVC model instrumented for surface 
temperature measurement using both thermocouples 
and thermal imaging technology. 

2) Record surface temperatures of the model in a subsonic 
wind tunnel environment. 

3) Compare the results from both measurement techniques in 
calculating the heat transfer coefficient, using 
Newton's Law of Cooling, and correlating the results 
with the theoretical solution for a flat plate under 
identical flow conditions. 

4) Utilize the thermal imaging technique to analyze the 
surface temperature map of an actual TVC vane 
configuration . 

5) Present conclusions and recommendations for future 
applications in this area. 
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II. MODEL DESIGN AND CONSTRUCTION 



A. MODEL DESIGN 

A double-wedge design was chosen for this experiment for 
several reasons, including ease of construction, similarity 
to typical TVC configurations used in previous research, and 
sufficient similarity to a flat plate to allow interpretation 
of results. A drawing of the wedge is included as Figure 1. 

The dimensions of the model were chosen to reflect essen- 
tially full-scale, but constrained by the necessity of having 
a Biot number less than 0.1. This constraint is mandated by 
the use of the lumped-heat capacitance analogy in this experi- 
ment to determine the heat transfer coefficient [Ref. 1] . 

This method of analysis assumes that the internal resistance 
of the body is negligible compared to the external resistance 
and, therefore, a uniform temperature distribution exists 
throughout the model. By ensuring that the model has a Biot 
number less than 0.1, this analytical method will produce 
reliable results. 

The convection heat loss from the heated model is assumed 
to equal the decrease in its internal energy. Thus, the 
energy balance expressing the rate of heat transfer, q, is: 

q = hA (T -TJ = “ c pV || 



where : 



h = heat transfer coefficient, BTU/hr-ft -°F 

2 

A = surface area, ft 

T = model temperature at time t, °F 

T^ = free stream temperature, °F 

c = heat capacity of model material, BTU/lbm-°F 

3 

V = model volume, ft 

k 5 thermal conductivity, BTU/hr-ft-°F 

3 

p = density, lbm/ft . 

Given the initial condition of T = T at time t = 0, the 

o 

solution to the differential equation becomes: 



T -T w 

^ — r^r - = exp - (Bi xFo) 

O ro 

where : 

Bi Biot number = 

Fo Fourier number 

B. CONSTRUCTION 

Copper was chosen as the model material due to its excel- 
lent thermal conductivity (k = 230 BTU/hr-ft- °F ) and reasonably 
good machinability . With the dimensions as shown in Figure 1, 

giving V/A = 0.01297, the surface heat transfer coefficient 

2 

could theoretically rise to well over h = 1400 BTU/hr-ft -°F 



hV 

kA 



= A(A) 2 t 

pc^V ; 
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before the restriction on Biot number would be violated. As 
the expected heat transfer coefficient was less than 60, these 
design decisions ensured excellent accuracy in applying the 
method . 

The copper wedge was instrumented with copper-constantin 
(type T) thermocouples, 30 gauge. The thermocouples were 
placed as shown in Figure 1, with one of the thermocouples 
embedded at depth of 0.30 inches to detect any distinguishable 
temperature difference between the surface and the interior. 

As can be seen in the subsequent tables and graphs, the assump- 
tion of Newtonian cooling proved to be accurate. The wedge 
was also covered with a thin coat of flat black paint to in- 
crease the emissivity of the surface to a minimum of 0.95 
[Ref. 1] in order to enhance the thermal image detection of the 
infrared system. 

The thermocouples were attached to a dual-channel strip 
chart recorder, and to a digital temperature sensor. The 
digital readout was used to obtain reference values of wedge 
temperature during heatup, and free stream temperatures during 
testing. All instruments and thermocouples were calibrated 
from 32 °F to 150 °F using the Department of Mechanical Engineer- 
ing calibration facility. 

The wedge was then mounted on two aluminum stanchions, 
attached to the plexiglass insert in the floor of the wind 
tunnel test section. Within the plexiglass, a six-inch diam- 
eter hole was cut for insertion of the infrared-transparent 
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window through which the thermal imaging system could view 
the wedge surface. This arrangement is shown in Figure 2. 
Plexiglass sidewalls were added to form a thermal barrier to 
attempt to negate any spanwise heat transfer from the wedge 
sides, and to form a channel for uniform air flow past the 
wedge. The model was heated by using an electric pad, capable 
of raising temperatures from ambient to approximately 140°F 
within a thirty minute time frame. 
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III. THERMAL IMAGING SYSTEM 



A. DESCRIPTION 

The thermal imaging system used in this procedure is a 
Probeye Thermal Video System , Series 4300, manufactured by 
Hughes Aircraft Company [Ref. 2], This system provides real- 
time temperature maps of objects through 16 color bands dis- 
played on a color monitor. The temperature range of the system 
is from -4°F to 536°F (-40°C to +280°C). Each color on the 
scale can represent temperature sensitivities from 0.5 to 
36°F (0.5 to 20°C) , depending upon the desired overall tempera- 
ture range. Emissivity adjustments can be made from 0.0 to 

1.0 in increments of 0.01, in the spectral wavelength band of 

2.0 to 5.6 ymeters. The unit includes a movable cursor which 
can be placed at any location on the viewed surface, giving a 
digital readout of the temperature at that location, and a 
real time clock/calendar precise to 0.1 seconds for determining 
time-variant conditions. 

The use of this imaging system also required the use of 
an infrared optic window for viewing the model while maintain- 
ing tunnel integrity. A six-inch diameter, two millimeter 
thick window fabricated from magnesium fluoride was obtained 
from Eastman Kodak Company. This type of material, tradename 
IRTRAN 1, was chosen to match the optical wavelength charac- 
teristics of the imaging system. As manufactured, the window 
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allowed over 95% transmittance in the 2.0 to 5.8 ymeter 
range, while capable of withstanding pressure differentials 
of 1 atmosphere at 25°C, with a factor of safety of four 
[Ref. 3] . 
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IV. EXPERIMENTAL PROCEDURE 



A. WIND TUNNEL TESTING 

The model assembly, as shown in Figure 2, was installed 
in the test section of the subsonic wind tunnel. This wind 
tunnel is a pull-down system capable of producing test section 
free stream velocities of up to 335 feet/second by means of a 
variable speed fan located at the tunnel exit. The free-stream 
turbulence level of the tunnel (0.2%-0.5%) was deemed suita- 
ble for the purposes of this experiment [Ref. 4] . 

At the start of each test, with the tunnel fan off, the 
wedge was wrapped with the heating pad, and brought up to 
temperatures of approximately 130 °F. The pad was then removed, 
wind tunnel access doors closed, and the fan started. Runs 
were made at 50% and 75% fan speed, corresponding to 200 and 
306 feet/second, respectively.. Free stream velocities were 
calculated using the installed manometer, calibrated in inches 
H^O. Free stream temperatures were determined by both a 
calibrated thermometer and a thermocouple inserted in the flow. 

The cooling process was then recorded using the installed 
thermocouple output to the strip chart recorders, and the ther- 
mal imaging system connected to a video tape recorder for 
simultaneous data acquisition. In initial tests, the tempera- 
ture sensitivity settings on the imaging system were selected 
at 5.0°F, which yielded an 80 degree range from 50°F to 130°F, 
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to view the entire cooling process. In subsequent tests, a 
sensitivity of 2.0°F was selected, giving a 32 degree range 
from 90°F to 122 °F. This lower sensitivity/range setting was 
necessary to obtain a good agreement with the thermocouple 
data . 

After a considerable amount of system testing, this process 
was undertaken four times at each of the two fan speeds for 
data evaluation and correlation. 

B. DATA REDUCTION 

Applying the formula for the lumped-capacitance analysis, 
the transient data from the thermocouple readings and the 
thermal imaging system were used to deduce the corresponding 
heat transfer coefficients. The calculations were performed 
using TKI SOLVER , a microcomputer spreadsheet [Ref. 5]. Tables 
1 and 2 exhibit the program equations and variables, which 
provide not only the heat transfer coefficient, but also an 
estimate of the uncertainty of that calculation, and the 
corresponding Biot number, for each of the two sensor types. 

The uncertainty analysis of the heat transfer coefficient 
calculation is included as Appendix A. 

The reduction of the thermocouple data was fairly straight- 
forward, as the output was temperature versus time. However, 
reducing the thermal imaging system data required running the 
video tape recording frame by frame (each frame corresponding 
to 1/60 second) to determine the temperature at a given time. 
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This method generated the largest uncertainty in the calculation 
because the visual color bands and digital readout of the sys- 
tem changed only when the temperature of the body changed 
more than the sensitivity setting. For example, at a 5.0°F 
setting, the readout would persist for perhaps three seconds 
of test time before changing from 110 °F to 105 °F, while in 
reality, as displayed by the thermocouple readings, the surface 
temperature could be any temperature within this range at a 
particular instant. (This problem is analogous to that due 
to a large interval size in a finite-difference calculation.) 
There also existed occasions when, particularly at temperatures 
near thermal equilibrium, the readout would take relatively 
long periods of time to stabilize, which made determining the 
exact temperature/time data difficult. However, these diffi- 
culties were overcome by operator experience and use of the 
narrower sensitivity setting for correlating data in smaller 
segments rather than over the entire cooling span. 
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V. EXPERIMENTAL RESULTS 



A. FLAT PLATE CORRELATION 

The analysis of a flat plate of identical material, length, 
and width in identical flow conditions was used to attain a 
degree of confidence in the accuracy of the results of the 
wedge experiment. This analysis is based upon determining 
the Reynolds and Prandtl numbers for the flow, and solving the 
following equation for the average Nusselt number: [Ref. 6] 

Nu = 0.664 Re 1//2 Pr 1//3 



The Nusselt number was then used to determine the average 
heat transfer coefficient from: 

Nu k 

h = £ 



where 



£ = length of flat plate, ft. 

These procedures resulted in a heat transfer coefficient of 
2 

22 BTU/hr-ft -°F for free stream velocity of 200 ft/sec, and 
2 

28 BTU/hr-ft -°F at velocity of 306 ft/sec. The results of 
this calculation are shown in Appendix B. 
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B. THERMOCOUPLE RESULTS 



Using the strip chart recorder traces, the heat transfer 

coefficient at the surface thermocouple location was determined 

from the lumped-capacitance solution of exponential temperature 

decrease with time. For the free stream velocity of 200 ft/ 

2 

sec, h varied from 27.5 to 33.2 BTU/hr-ft -°F, with an average 

2 

uncertainty band of ±1.75 BTU/hr-ft -°F. At velocity of 

306 ft/sec, h had values ranging from 41.03 to 42.7 BTU/ 

2 2 
hr-ft -°F, with average . uncertainty of ±1.80 BTU/hr-ft -°F. 

These calculations are shown in Tables 4-11. 

C. THERMAL IMAGING SYSTEM RESULTS 

As shown in Tables 12-19, using the same analysis and 

temperature range, h values were determined to be 27.7 to 

33.5 BTU/hr-f t^-°F, with an uncertainty band of ±1.9 BTU/ 

2 

hr-ft -°F for the lower free stream velocity, while values of 
39.9 to 41.8 were obtained at higher speed. Figures 3 and 4 
show the comparison of wedge surface temperature versus time 
for the two sensors at each velocity, and Figures 5 and 6 
show the calculated values from each sensor versus time, again 
at each speed. 

It is important to note that this measurement corresponds 
to conditions at a particular point on the wedge, that is, the 
location of the corresponding thermocouple. Further testing, 
with the cursor positioned at other points, would allow the 
development of a complete heat transfer map of the wedge sur- 
face. However, viewing of the entire wedge surface showed that 
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the transient isotherm variation was essentially uniform over 
a large portion of the surface (with the exception of narrow 
bands in the three-dimensional boundary layer regions near the 
sidewalls) . Figure 7 is a still photograph of the wedge taken 
from the video output, which clearly exhibits this condition. 
Thus it may be deduced that the values determined above, at 
the thermocouple - location , are representative of conditions 
over a major portion of the wedge surface. The ability to make 
such observations is a significant virtue of the infrared 
thermography method. 

D. TVC VANE ANALYSIS 

An actual TVC vane was received from NWC China Lake, Ca 
for analyzing the effect of complex geometry on convective heat 
transfer using the thermal imaging system. 

This vane, shown in Figure 8, was constructed of copper- 
impregnated tungsten (80% W, 20% Cu) , with the following 
material properties: 



= 1076.544 



lbm 

ft 3 



= 148 



BTU 

hr-f t-°F 



0.037 



BTU 

lbm-°F 



-5 3 

The vane volume was calculated to be 2.295 x 10 ft , with 

-1 2 

a corresponding surface area of 1.5551 x 10 ft . All values 
were then entered into a formula sheet as before, shown in 
Table 3. 

The vane was mounted on two aluminum stanchions with the 
same baseplate used by China Lake in their experiments. The 
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thermal imaging system was set up to view the vane through the 
side port of the subsonic wind tunnel with the IRTRAN optic 
window installed. The identical heating procedure was used, 
however, for reasons unknown, the pad would only reach 105 °F. 
This was deemed sufficient, as it still gave a 40 degree tem- 
perature differential. The wind tunnel was operated at 50% 
flow speed (200 ft/sec) . 

Three test runs were performed, with the cursor at the 
midspan for runs 1 and 3, and on the leading edge at the base- 
plate boundary for run 2. The cursor was moved for run 2 to 
investigate and analyze a more rapid temperature gradient 

noticed in that area in the previous run. Run 1 yielded an 

2 

average h of 41.5 BTU/hr-ft -°F ±3.23, run 2 yielded an average 
of 54.8 ±4.41, and run 3 gave an average of 42.4 ±3.9. The 
largest Biot number was calculated as 0.0063, well within the 
lumped-capacitance analysis criteria. 

Figure 9 is a color thermograph of the vane during cool- 
down. As compared to the thermogram of the double wedge 
(Figure 7), the temperature map is fairly uniform over the 
body, with the exception of faster cooling at the leading and 
trailing edges, and an interesting area near the top of the 
vane below the apex. While vane thickness would explain the 
shape of most of the isotherms, the upper region undoubtedly 
is affected by a unique flow field, which increases the heat 
transfer coefficient by approximately 8% (Table 22). Tables 
20-22 show the results of these tests, and Figures 10 and 11 
are graphs of temperature and h versus time, respectively. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The results obtained from the thermal imaging system are, 
within a small experimental uncertainty, in agreement with 
the results of the thermocouples for all conditions tested 
here. Both experimental method's results were also suffi- 
ciently close to the theoretical flat plate solution to allow 
confidence in their accuracy, primarily due to the small wedge 
angle used. 

The primary disadvantage of the thermal imaging system is 
the increase -in size of the uncertainty band as the tempera- 
ture sensitivity is widened. For large sensitivities of 5 
degrees or more, the results become comparatively less accurate, 
again similar to using large step sizes in finite element 
calculations. Future users of this system in experiments 
similar to this one will have to adjust the sensitivity settings 
with respect to the overall temperature range desired, most 
likely by trial and error, to get the degree of accuracy 
required . 

Nevertheless, the results obtained in this experiment 
lead to the conclusion that the thermal imaging system, when 
properly employed, would be a very good alternative to thermo- 
couples in applications similar to this one. It provides 
accurate local and global readings with no need for 
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complicated data acquisition/reduction equipment and extensive 
instrumentation. The system is transportable, very reliable, 
and easy to operate. Although the system requires special 
infrared optic material for viewing models in an enclosed 
environment, several vendors were found who are capable of 
providing the material in various shapes and sizes at a reason- 
able expense in time and money. 



B. RECOMMENDATIONS 

The findings of this research lead to a recommendation for 
using this system for surface temperature measurements in the 
next phase of this project, using the wedge model in a super- 
sonic wind tunnel for the same heat transfer coefficient calcu- 
lation. This model was designed for ease of installation in 
the supersonic tunnel, when it becomes available, and an optic 
window has been obtained that will allow infrared viewing of 
the model while withstanding the environment (temperature and 
pressure differentials) in supersonic flow. Particular problems 
that might be encountered in supersonic applications include: 

1) Sufficient difference in model and free stream 
temperatures to cause adequate heat transfer — the 
model may have to be heated prior to testing by 
some internal means (installation of wires for 
Joulean heating) to increase the recovery factor. 

2) The on-site supersonic tunnel has a significantly 
smaller run-time available due to its pressure 
vessel configuration, thus yielding short run 
durations (a few minutes) and long (hours vice 
minutes) pressurization durations between tests. 

3) As other model shapes and materials are considered, 
the necessity for maintaining the Biot number less 
than 0.1 could become critical, or another method 
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of analysis would have to be chosen if the lumped- 
capacitance assumption is not valid. 

It is also recommended that the next experimenter familiar- 
ize himself with the system thoroughly, and invest time in 
viewing the video tape outputs of the system to get experience 
in temperature determination. This experience will prove to 
be a great factor in enhancing the accuracy of the results. 
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APPENDIX A 



UNCERTAINTY ANALYSIS CALCULATIONS 



T -T . 

00 hAt \ 

0 = ^ = exp - (— ^) Formula used in calculation 

O 00 ^ 

, - J£n0(p c V) -1 £ . 

h = A~t ' solving for h 

h = F(0,t,A,V,p,c) 



au _ r/3h . n > 2 , . 3h . . 2 , . 3h . _ , 2 , ,3h .,.,2 

Ah " [( 30 A0) + { 3t At) + ( 3A AA) + 3V AV 



, 3h . ,2 , . 3h . .2. 
+ lTp" Ap) + 3c" AC) J 



1/2 



where A's are the uncertainty of the measured quantities. 



CL) || A0 = 



pcV r , N .2 , .1 ..,>2, 

£_[(. AD ) + AN) ] 



1/2 



where : 



2 2 

N = T - T , D = T -T , AN = (AT + AT ) , 

00 ' o 00 00 ' 



2 2 1/2 

AD = (AT + AT ) , AT = 1.0 °F, AT 

o 00 00 o 



= 1.0 °F 



AT = 2.0 °F for IR data, AT = 1 . 0 °F for thermocouple 

data 
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(2) || At = 


Jin 0 pcV AJ _ 
At 



where : 



At = 


= 0.1 seconds = 2.7778 E -5 hours 


from IR system 


timer precision. 


< 3 > lx iA ■ 


. - in V c V AA 

A Z t 



where : 



AA = 


= 0.001 in 2 = 6 . 9444E-6f t 2 


from machining 


precision . 


(4) f^AV = 


■ - in A ®t pc av ' 



where : 



AV = 


= 1.0E-9 in 3 = 5 . 78 78E-13 ft 3 


from machining 


precision . 


, c . 3h . 

(S) ^ Ap . 


Jin 0c V . 

’ ' At AP ' 



where : 



Ap = 


= 6.9082 lb ™ [Ref. 7] 

ft 
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( 6 ) 



( 6 ) ^7 Ac 


£n 0 p V 

"'At A ' 



where : 



Ac 


- °- 001 lbm-°F !Ref ' 7J 


All above 


equations were included in the data reduction 



programs to solve for the uncertainty in h (delh) . 
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APPENDIX B 



FLAT PLATE SOLUTION 



Given: T^ = 60°F = 520°R 1"H 2 0 = 0.0361 psi 



Find: h from the empirical flat plate in parallel flow 

equation : 



Nu = 0.664 Re 1//2 Pr 1//3 



x 



[Ref. 6] 



h = 



Nu k 



(1) Re : For T = 520°R: [Ref. 7] 

X 00 



air 



= 1.2039E-5 



lbm 



f t-sec 



Pr = 0.71 k = 0.01561 



BTU 



hr-f t-°F 



y . = 0.0743 

1 3 1 y* 



air 



lbm 

ft 3 



x = 2.50 in = 0.2083 ft 



U 



2 _ Ap 2 g 

oo y 



(8.9"H 2 0) (0.0361) (64.348) (144) 
0T0743 



U 2 = 40080.4 ft 



ft 



U^ = 2 00.2 — for experimental 

low speed flow 



For Ap = 



20.8” H n O, U 

2 oo 



= 306.1 — for experimental 

high speed flow 
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YU x 



Re 



“ _ (0.0743) ( 200 . 2 ) (0.2083) 



1.2039E-5 



= 2.573E5 



Re 



(0 . 0743) (306.1) (0 . 2083) 



1. 2039E-5 



= 3.935E5 



Assuming Re cr ^ t = 5.0E5, laminar flow conditions exist in both 
cases . 



(2) Nu: For low speed flow: 



Nu = 0.664 (2.574E5 ) 1//2 (0.71) 1/3 ) = 300.5 



For high speed flow: 



Nu = 0.664 (3.935E5) 1//2 (0.71) 1/3 = 371.6 



(3) h: For low speed flow: 



- Nu k (300.5) (0.01561) 

h — ~ — 



0.2083 



= 22.5 



BTU 



hr-ft - °F 



For high speed flow: 



h = Nu k = (371.6) (0 .01561) 27 8 



x 



0.2083 



BTU 

hr-f t 2 -°F 
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APPENDIX C 



TABLES 



TABLE 1 

FORMULA SHEET FOR THERMOCOUPLE DATA REDUCTION 



St 


Input 


Name 


Output Unit 


Comment 


L 




theta 




dimensionless 

temperature 


L 




Tt 


F 


temperature at time t 






Tinf 


F 


free stream 
temperature 






Ti 


F 

btu 


initial wedge 
temperature 


T. 




h 


heat transfer 


XJ 




2 








hr-ft -°F 


coefficient 


T 




delh 

a 


btu 


uncertainty in h value 
wedge wetted perimeter 


±j 


.27777778 


hr-f t 2 -°F 
ft A 2 


L 




t 


hours 






558 


ro 


lbm/f t~ 3 


copper density 




.091 


c 


btu/lb-F 


specific heat capacity 




.00346013 


V 


ft~3 


wedge volume 


L 




q 


btu/hr 


heat transfer 




230 


k 


btu 

hr-f t-°F 


wedge conductivity 


L 




BIOT 






L 




time 


seconds 








u 




(dh/dtheta) * (delta theta 






vv 




(dh/dt) * (delta t) 






w 




(dh/dA) * (delta a) 






X 




(dh/dV) * (delta V) 






Y 




(dh/dro) * (delta ro) 






z 




(dh/dc )* (delta c) 
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TABLE 1 (CONTINUED) 



S Rule 



* theta 



(Tt-Tinf) 
(Ti-Tinf ) 



* theta 



exp 



(- (h*a*t ) 
(ro*c*V) ) 



* q = 

* BIOT 



* time 



h*a* (Tt-Tinf) 

h*V 

(k*a) 

= t*3600 



* u = (ro*c*V)/ (theta*at) )*( (1 .4142* (Tt-Tinf )/ (Ti-Tinf ) ~2) 

+ (1 . 8868/ (Ti-Tinf )) "2) 

* vv = (2. 777777 e-5*ln (theta ) *ro*c*V/ (a*t ~2 ) ) 

* w = ( In ( theta ) *ro*c*V* 6 . 9444e-6/ (a~2*t ) ) 

* x = (In (theta) *ro*c*5. 787e-13/(a*t) ) 

* y = ( In ( theta) *c*V*6.908 2/ (a*t) ) 

* z (In (theta) *V*ro* . 001/ (a*t) ) 

* delh = sort (u / '2+w / '2+w A 2+x / '2+y / '2+z /v 2) 
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TABLE 2 



FORMULA SHEET FOR INFRARED SYSTEM DATA REDUCTION 



St 


Input 


Name 


Output Unit 


Comment 


L 




theta 




dimensionless temperature 


L 




Tt 


F 


temperature at time t 






Tinf 


F 


free stream temperature 






Ti 


F 


initial wedge temperature 


L 




h 


btu 


heat transfer coeffi- 








hr-f t^-°F 


cient 


L 




delh 


btu 


uncertainty in h value 
wedge wetted perimeter 




.27777778 


a 


hr-f t4-°F 
ft 2 


L 




t 


hours 






558 


ro 


lbm/ft~3 


copper density 




.091 


c 


btu/lb-F 


specific heat capacity 




.00346013 


V 


ft~3 


wedge volume 


L 




q 


btu/hr 


heat transfer 


L 


230 


k 

BIOT 


btu 

hr-f t-°F 


wedge conductivity 


L 




time 


seconds 








u 




(dh/dtheta) * (delta theta) 






vv 




(dh/dt) * (delta t) 






w 




(dh/dA) * (delta a) 






X 




(dh/dV) * (delta V) 






y 




(dh/dro )* (delta ro) 






z 




(dh/dc) * (delta c) 



S Rule 



* theta 



(Tt-Tinf ) 
(Ti-Tinf ) 



* 



theta = exp ( 



- (h*a*t) , 
(ro*c*V) 



* q = h*a* (Tt-Tinf) 
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TABLE 2 (CONTINUED) 



* BIOT 

* time 

* u = 

* vv = 

* w = 

* x = 

* Y = 

* z = 

* delh 



h*V 

(k*a) 

= t *3600 

(ro*c*V/ (theta*a*t) ) * ( (1 . 4142* (Tt-Tinf )/ (Ti-Tinf ) "2 ) 
+ (2 . 2361/ (Ti-Tinf ) ) ~2) 

(2. 777777 e- 5* In (theta) *ro*c*V/ (a*t~2 ) ) 

(In (theta) *ro*c*V*6 .9444e-6/ (a~2*t) ) 

(In (theta) *ro*c*5 . 787e-13/ (a*t) ) 

(In (theta) *c*V*6 . 9082/ (a* t) ) 

(In (theta) *V*ro* . 001/ (a*t) ) 

= sqrt (u'' s 2+vv / '2+w / '2+x'' v 2+y'''2 + z / '2) 
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TABLE 3 



FORMULA SHEET FOR TVC VANE DATA REDUCTION 



S Rule 
* theta 



(Tt-Tinf ) 
(Ti-Tinf) 



* 



theta = exp ( 



- (h*a*t) > 
(ro*c*V) 



* 



* 



q = 

BIOT 

time 



h*a* (Tt-Tinf) 

h*V 

(k*a) 

= t*3600 



* u = (ro*c*V/(theta*a*t) ) * ( (1.4142* (Tt-Tinf) / (Ti-Tinf) ~2) 

+ (2 . 2361/ (Ti-Tinf ) ) ^2) 

* vv = (2 . 77 7 7 77 e- 5* In’ (theta) *ro*c*V/ (a*t~2) ) 

* w = (In (theta) *ro*c*V*6 . 9444e-6/ (a~2*t) ) 

* x = (In (theta) *ro*c*5 . 787e-13/ (a*t) ) 

* y = (In (theta) *c*V*6 . 9082/ (a*t) ) 

* z = ( In ( theta) *V*ro* . 001/ (a*t) ) 

* delh = sqrt (u''2+w / '2+w / '2+x / '2+y /v 2+z , '2) 



35 



TABLE 3 (CONTINUED) 



St 


Input 


Name 


Output Unit 


Comment 


L 




theta 




dimensionless temperature 


L 


81 


Tt 


F 


temperature at time t 




61 


Tinf 


F 


free stream temperature 




97 


Ti 


F 


initial TVC vane 


L 




h 


btu 


temperature 

heat transfer coefficient 






L 




delh 


hr-ft - °F 


uncertainty in h value 




.155507 


a 


ft" 2 


TVC vane surface area 


L 




t 


hours 






1076.544 


ro 


lbm/f t~ 3 


TVC vane density-80% W, 
20% Cu 




.037 


c 


btu/lb-F 


specific heat capacity 




.00229545 


V 


ft A 3 


TVC vane volume 


L 




q 


btu/hr 


heat transferred 




148 


k 


btu 

hr-ft-°F 


TVC vane conductivity 


L 




BIOT 






L 


47.3 


time 


seconds 
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TABLE 4 



RESULTS OF THERMOCOUPLE CALCULATIONS, LOW TUNNEL SPEED 



Run : 


#1 


Tunnel Air Velocity: 


: 201.3 ft 

sec 


Average Biot Number : 


: 0.00164 


Average Uncertainty 


in h: ±2.28 

hr-ft -°F 



Time (Seconds) 


Temperature (°F) h tj 1 r-ft-°F^ 


0.0 


117.5 


6.7 


112.5 30.92 


14.5 


107.5 30.00 


22.7 


102.5 30.32 


32.0 


97.5 30.42 


43.0 


92.5 30.21 


55.0 


87.5 30.54 
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TABLE 5 



RESULTS OF THERMOCOUPLE CALCULATIONS, LOW TUNNEL SPEED 



Run : 


#2 


Tunnel Air Velocity: 


: 193.3 

sec 


Average Biot Number: 


: 0.00156 


Average Uncertainty 


in h: ±1.95 bt ^ 

hr-ft -°F 



Time (Seconds) 


Temperature (°F) h ^r-ft- 0 F^ 


0.0 


114 . 5 


3.5 


112.5 24.32 


5.6 


110.5 30.99 


9.3 


108.5 28.56 


12.8 


106.5 28.24 


16.0 


104.5 28.85 


21.0 


102.5 26.97 


24.6 


100.5 27.48 
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TABLE 6 



RESULTS OF THERMOCOUPLE CALCULATION, LOW TUNNEL SPEED 



Run : 


#3 


Tunnel Air Velocity: 


: 197.6 ft 

sec 


Average Biot Number: 


: 0.00159 


Average Uncertainty 


in h: ±1.53 bt ^ 

hr-ft -°F 


Time (Seconds) 


Temperature (°F) h (fckft °F^ 


0.0 


120.0 


10.0 


113.3 28.70 


20.0 


107.0 29.47 


30.0 


101.7 29.39 


40.0 


96.7 ' 29.92 


50.0 


92.8 29.54 


60.0 


89.2 29.54 
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TABLE 7 



RESULTS OF THERMOCOUPLE CALCULATIONS, LOW TUNNEL SPEED 



Run : 

Tunnel Air Velocity: 
Average Biot Number: 
Average Uncertainty in h: 



#4 

207.9 -Q- 
sec 

0.00182 



± 2.41 



btu 

hr-f t 2 - °F 



Time (Seconds) 

0.0 

5.0 

10.0 

15.0 

20.0 

25.0 

30.0 

35.0 

37.0 



Temperature (°F) 

111.0 

107.7 

104.3 

101.4 

98.7 
9 6.2 

93.7 
91.6 
91.0 



h ( hr-ft-°F ) 



32.44 

34.23 

33.72 

33.71 

33.69 

33.82 

33.57 

33.17 
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TABLE 8 



RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 



Run : 




#5 






Tunnel Air Velocity: 


308.2 


ft 

sec 




Average 


Biot Number: 


0.0022 






Average 


Uncertainty 


in h : ±1.88 


btu 




hr-f t 2 -°F 




Time (Seconds) 


Temperature ( 


°F) h 


, btu 
hr-f t- °: 


0.0 




119.0 






5.0 




114.3 




40.37 


10.0 




109.7 




41.45 


15.0 




105.5 




41.73 


20.0 




101.7 




41.95 


25.0 




98.3 




42.03 


30.0 




95.2 




42.00 


35.0 




92.2 




42.26 


37.0 




91.0 




42 . 66 
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TABLE 9 



RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 



Run : 


#6 






Tunnel Air Velocity: 303.8 


ft 

sec 




Average 


Biot Number: 0.0022 






Average 


Uncertainty in h: ±2.11 


btu 




hr-ft 2 -°F 




Time (Seconds) Temperature ( 


°F) h 


, btu 
hr-f t- 


0.0 


119.0 






5.0 


114.5 




38.15 


10.0 


110.3 




38.64 


15.0 


106.7 




37.48 


20.0 


102.1 




40.84 


25.0 


98.5 




41.52 


30.0 


95.6 




41 . 07 


35.0 


92.4 




41.88 


38.0 


91.0 




41.54 
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TABLE 10 



RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 



Run : # 7 



Tunnel Air Velocity: 


306.0 ft 
sec 




Average Biot Number: 


0.0022 




Average Uncertainty 


in h: ± 2.28 

hr-ft -°F 




Time (Seconds) 


Temperature (°F) h 


, btu 
l hr-ft- 


0.0 


119.0 




5.0 


114.4 


38.06 


10.0 


110.0 


39.13 


15.0 


105.7 


40.30 


20.0 


102.0 


40.32 


25.0 


98.5 


40.70 


30.0 


95.4 


40.57 


35.0 


91 . 4 


43.12 


37.5 


91.0 


41 . 03 
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TABLE 11 



RESULTS OF THERMOCOUPLE CALCULATIONS, HIGH TUNNEL SPEED 



Run : 


#8 






Tunnel Air Velocity: 


306.0 


ft 

sec 




Average Biot Number: 


0.0022 






Average Uncertainty 


in h : ±1.96 


btu 




hr-f t 2 -°F 




Time (Seconds) 


Temperature ( 


°F) h 


, btu 

hr-f t- 


0.0 


119.0 






5.0 


114.4 




39.57 


O 

o 
1 — 1 


110.0 




39.89 


15.0 


105 . 8 




40.81 


20.0 


102.2 




40.61 


25.0 


98.6 




41.18 


30.0 


95.5 




41.23 


35.0 


92.5 




41.70 


38.0 


91.0 




41.53 
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TABLE 12 



RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW TUNNEL SPEED 



Run : 

Tunnel Air Velocity: 
Average Biot Number: 
Average Uncertainty in h: 



#1 



201.3 



ft 

sec 



0.00164 



± 1.88 



btu 

hr-f t 2 -°F 



Time (Seconds) 



Temperature (°F) 



h ( 



btu 

hr-f t- °F 



) 



0.0 

6.1 

13.0 

22.0 

31.5 
41.7 

52.5 



117.5 




112.5 


33.96 


107.5 


33.46 


102 . 5 


31.29 


97.5 


30.90 


92.5 


31.15 


87 . 5 


31.99 
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TABLE 13 



RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW TUNNEL SPEED 
Run : 

Tunnel Air Velocity: 

Average Biot Number : 

Average Uncertainty in h: 

j— ^ 

Time (Seconds) Temperature (°F) h • ( fc ~ r lf ^_ ~ o ~p) 



#2 



193.3 



ft 

sec 



0.00151 



± 2.16 



btu 



hr-ft - °F 



0.0 


114.5 




3.1 


112.5 


28 .40 


6.1 


110 . 5 


28 .46 


9 .5 


108 .5 


27 .96 


13 .3 


106 . 5 


27 .18 


17.0 


104 .5 


27 .15 


21 .5 


102 . 5 


27 .33 


24.4 


100.5 


27 . 71 
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TABLE 14 



RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW 



Run : 

Tunnel Air Velocity: 
Average Biot Number: 
Average Uncertainty in h: 



#3 



197.6 



ft 

sec 



0.00161 



± 1.98 



btu 



hr-ft -°F 



Time (Seconds) Temperature (°F) 



0.0 


119.0 


5.5 


115.0 


11.5 


111.0 


18.9 


107.0 


26.5 


103.0 


32.3 


99.0 


40.4 


95.0 


50.4 


91.0 



TUNNEL SPEED 



(— ) 

v hr-f t- °F ; 



30.68 

30.52 

29.05 

28.91 

31.15 

31.54 

31.32 
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TABLE 15 

RESULTS OF INFRARED SYSTEM CALCULATIONS, LOW 



Run : 

Tunnel Air Velocity: 
Average Biot Number: 
Average Uncertainty in h: 



#4 

207.9 -££- 
sec 

0.0018 



± 2.14 



btu 

hr-f t 2 -°F 



Time (Seconds) 

0.0 

2.8 

6.0 

9.4 

12.4 

16.3 

19.3 

22.8 

27.0 
31.8 
36.6 



Temperature (°F) 

111.0 

109.0 

107.0 

105.0 

103.0 

101.0 

99.0 

97.0 

95.0 

93.0 

91.0 



TUNNEL SPEED 



, btu 
^hr-f t- °F 



34 . 61 
33.02 
32.35 
33.48 

32.64 
33.94 
34.44 
34.19 

33.65 
33.53 
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TABLE 16 



RESULTS OF INFRARED SYSTEM CALCULATIONS, HIGH TUNNEL SPEED 



Run : 

Tunnel Air Velocity: 
Average Biot Number : 
Average Uncertainty in h: 



#5 



308.2 



ft 

sec 



0.0021 



± 1.85 



btu 

hr-f t 2 -°F 



Time (Seconds) Temperature (°F) 



h ( 



btu 

hr-f t-°F 



0.0 


119.0 


4.7 


115.0 


8.9 


111.0 


13.7 


107.0 


19.3 


103.0 


23.8 


99.0 


29.6 


95.0 


37 . 8 


91.0 



35.91 

39.44 

40.08 

39.70 

42.27 

43.05 

41.75 
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TABLE 17 



RESULTS OF INFRARED SYSTEM CALCULATIONS, HIGH TUNNEL SPEED 



Run : 

Tunnel Air Velocity: 
Average Biot Number: 
Average Uncertainty in h: 



#6 

303.8 

sec 

0.0018 



± 1.77 



btu 

hr-ft 2 -°F 



Time (Seconds) 



Temperatur e ( °F) 



h ( 



btu 

hr-f t- °F 



) 



0.0 


119 . 0 


5.9 


115.0 


10.8 


111.0 


16.2 


107.0 


21.6 


103.0 


26.5 


99.0 


32.7 


95.0 


39 . 6 


91.0 



35.45 
37 . 69 
36.91 
35.47 

37.97 

38.97 
39.86 
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TABLE 18 



RESULTS OF INFRARED SYSTEM CALCULATIONS, HIGH TUNNEL SPEED 



Run : 



#7 



Tunnel Air Velocity: 



306.0 



ft 

sec 



Average Biot Number: 



0.0021 



Average Uncertainty in h: ±1.86 



btu 



hr-ft -°F 



Time (Seconds) 



Temperature (°F) 



h ( 



btu 

hr-f t- °F 



) 



0.0 


119.0 


5.0 


115.0 


9.0 


111.0 


14.0 


107.0 


18.9 


103.0 


24.0 


99.0 


29.7 


95.0 


36.8 


91.0 



36.14 
38.26 
38.45 
39 . 70 
41.00 
41.90 
41.81 
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TABLE 19 



RESULTS OF INFRARED SYSTEM CALCULATIONS, HIGH 
Run : 

Tunnel Air Velocity: 

Average Biot Number : 

Average Uncertainty in h: 

Time (Seconds) Temperature (°F) 



#8 



306.0 



ft 

sec 



0.0022 



± 1.83 



btu 



hr-ft -°F 



0.0 


119.0 


4.9 


115.0 


10.3 


111.0 


14.7 


107.0 


20.6 


103.0 


25.9 


99.0 


31.5 


95.0 


38.2 


91.0 



TUNNEL SPEED 



(— k tu } 

'hr-f t- °F^ 



37.14 

37.46 

37.45 
37.20 
38.85 

40.45 
41.32 
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TABLE 20 



RESULTS OF TVC VANE CALCULATIONS 



Run : 



Tunnel Air Velocity: 



Average Biot Number: 



Average Uncertainty in h: 



#1 



214.3 



ft 

sec 



0.0041 



± 3.23 



btu 



hr-ft - °F 



Time (Seconds) Temperature (°F) 



h ( 



btu 

hr-f t-°F 



) 



0.0 


99.0 


5.5 


95.0 


10.8 


91.0 


17.0 


87.0 


23.5 


83.0 


31.6 


79.0 


42.2 


75.0 



38.52 

41.41 

41.89 

43.19 

43.31 

42.50 



53 



TABLE 21 



RESULTS OF TVC VANE CALCULATIONS 



Run : # 2 

Tunnel Air Velocity: 220.5 

1 sec 

Average Biot Number: 

Average Uncertainty in h: 



Time (Seconds) Temperature (°F) h tL o p ) 



0.0 


97.0 




4.0 


93.0 


62.33 


9.0 


89.0 


59.10 


16.4 


85.0 


52.33 


23.2 


81.0 


53.63 


33.0 


77.0 


52.01 


47.3 


73.0 


49.17 



0.0053 



± 4.41 



btu 



hr-ft -°F 
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TABLE 22 



RESULTS OF TVC VANE CALCULATIONS 



Run : 


#3 


Tunnel Air Velocity: 


: 220.5 ft 

sec 


Average Biot Number : 


: 0.0042 


Average Uncertainty 


in h: ±3.90 bt ^ 

hr-ft -°F 


Time (Seconds) 


Temperature (°F) h (fokft- °F^ 


0 . 0 " 


97.0 


6.2 


93.0 41.43 


13.0 


89.0 42.25 


21.7 


85.0 40.95 


28.7 


81.0 45.06 


40.4 


77.0 44.39 
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APPENDIX D 



FIGURES 
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Figure 1 . Diagram of Wedge Model 
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Figure 2. Rear View of Wedge Model Assembly 
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TEMP VS TIME-LOW SPEED 
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Figure 3. Graph of Wedge Temperature vs Time, Low Speed Flow 



TEMP VS TIME-HIGH SPEED 
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Figure 4. Graph of Wedge Temperature vs Time, High Speed Flow 
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Figure 5. Graph of Wedge h vs Time, Low Speed Flow 
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Figure 7. 



Color Thermogram of Wedge 
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Figure 8. Photograph of TVC Vane 





Flow Direction 
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TVC VANE TEMPERATURE VS TIME 
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Figure 10. Graph of TVC Vane Temperature vs Time 
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Figure 11. Graph of TVC Vane h vs Time 
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